Numerical study of a superconductor-insulator transition in a half-filled Hubbard 

chain with distant transfers 
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The ground state of a one-dimensional Hubbard model having the next-nearest neighbor hopping 
{t') as well as the nearest-neighbor one (f) is numerically investigated at half-filling. A quantum 
Monte Carlo result shows a slowly decaying pairing correlation for a sizeable interaction strength 
([/ < 2t), while the system is shown to become insulating for yet larger U > Uc ~ 3t from a 
direct evaluation of the charge gap with the density-matrix renormalization group method. The 
results are consistent with Fabrizio's recent weak-coupling theory which suggests a transition from 
a superconductor into an insulator at a finite U . 
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For the past several years, the existence of a gap in 
the spin excitation has been suspected to be a key in 
the superconductivity in one-dimensional (ID) strongly 
correlated systems. This vieW|-|lijis motivated inten- 
sive studies on ladder ajystemsjj13 and systems with 
dimerizedlj or frustratedu spin- spin interactions. Some 
of these studies have revealed that the existence of a 
spin gap can indeed lead to a dominance of pairing 
correlations even in purely repulsive models. More re- 
cently, Fabrizio studied a ID Hubbard model having the 
next-nearest neighbor (NNN) hoppings in addition to 
nearest- neighbor (NN) ones. Such a model may have 
some relevance to the ID cuprates with large ratios of 
NNN to NN transfers, such as SrCu02 having a zigzag 
structureji^ or the roecntly discovered superconductor, 
Sro.4Cai3.6Cu2404i.84ll3, where layers of 2-leg ladders al- 
ternate with layers of ID chains with 90° Cu-O-Cu bond 
angles. Making use of a weak-coupling renormalization, 
Fabrizio has predicted the existence of a spin gap in a 
certain parameter regime and also a 'superconductor- 
insulator transition' (or more precisely a transition from 
a metal with dominating pairing correlation to an insu- 
lator since we are talking about ID systems) suLa certain 
Hubbard interaction U = Uc at half-fiUingJlj A tran- 
sition from a superconductor directly into an insulator 
is intriguing, but it is beyond the scope of the weak- 
coupling theory to establish where (i.e., for which finite 
value of the interaction strength) the transition actually 
takes place. 

The purpose of the present study is to look into 
the transition numerically, and to actually evaluate Uc- 
Our strategy here is to employ two complementary ap- 
proaches: the quantum Monte Carlo (QMC) method for 
small t/'s and the density-matrix renormalization group 
(DMRG) for large C/'s. Combination of the two has in- 
deed enabled us to establish the existence of such a tran- 
sition at Uc ^ it. 

The Hamiltonian under consideration is given, in stan- 
dard notations, as 



^ = -^ E(4c.+ia + h.c.) - t' ^(4c,+2. + h.c.) 
-fC/^n.tn^, (1) 



where t and t' are NN and NNN hoppings, respectively, 
and U is the Hubbard repulsion. We focus our attention 
to half- filling, n = 1, where the signs of t and t' are irrel- 
evant due to an electean-hole symmetry. Here we shall 
take t > and i' < O.B 

Let |-u£ first recapitulate Fabrizio's weak-coupling 
theory.E] If \t' /t\ is small enough, the system can be 
mapped to the t' — Q case, because there is no essen- 
tial difference in the vicinity of the Fermi level. For 
\t'\ > tc{n), on the other hand, the Fermi level inter- 
sects the one-electron band at four /c-points (±fc^ and 
±kp), where tc{n) is a function of the band filling (with 
tc(l) — O.bt). Namely, there are two right-moving and 
two left-moving branches in the terminology of the weak- 
coupling Tomonaga-Luttinger (TL) theory. When the 
Umklapp processes are absent, the situation is, as far as 
the weak-coupling picture is concerned, essentially iden- 
tical to the two-leg Hubbard ladder, where a spin gap 
opens when the two Fermi velocities are not too differ- 
ent. Thus a spin gap also opens in the present model for 
a certain range of the ratio of the Fermi velocities, i.e., for 
\t'\ > i^(n)(> tc{n)). In the presence of a spin gap, the 
pairing correlation dominates for small enough U ii 
two-leg ladder when doped (i.e., at non-half- filling)) 
and this is also the case with the present NNN model. 

An important difference, however, appears at half- 
filling, where Umklapp processes emerge. In the two-leg 
ladder, in which the Umklapp process is a two-electron 
scattering process at half-filling, an infinitesimal U is 
enough to make the Umklapp process relevant in the 
renormalization, leading to an opening of a charge gap. 
In contrast, Umklapp processes are only higher-order at 
half-filling in the present NNN model, in which a four- 
electron scattering is involved. 

The dominant pairing considered in Fabrizio's weak- 



coupling picture is 



Cfei T'^-fcU ~ cj,2 -fC_fe2 1, 



(2) 



which has the correlation function of the form ^ i/r^f^^ 
(r: the real space distance). This pairing dominates over 
a charge-density wave (or dimer wave in Fabrizio's termi- 
nology) whose correlation is ~ 1/r^^ . If we assume that 
the analysis may be extended to finite U^s, the renormal- 
ization equation has such a structure that K > 0.5 for 
U < Uc (with iiT ^ 1 for [/ — > -1-0), where the pair- 
ing correlation is the most slowly- decayiiig|-(~ l/r"-^ in 
the^ weak-coupling limit as in doped two-a □ and three- 
legj'Q Hubbard ladders). As U is increased the Umk- 
lapp process becomes relevant to make the system insu- 
lating precisely when the two functional forms coincide 
(l/ri/2^ = l/r2^, i.e., at [/ = Uc)- 

Now let us turn to our numerical analysis for finite [/'s. 
We take t' ~ — 0.8i as a typical value of t' > tc = 0.5t. 
We first present the QMC result of the pairing correlation 
function in the ground state. We have adopted the pro- 
jector Monte Carlo method, in which we implemeni—jtlie 
stabilization procedure adopted by several authors .EtEEI 
We have taken projecting time r of up to ~ 60/t with 
Trotter slices L of t/L < 0.2/U to assure the conver- 
gence of the correlation function. The details of the 
calculation is similar to that for the two-leg Hubbard 
ladder in ref. ^ We take 60 electrons on 60 sites with 
t' = —0.78t. For this set of parameters, the energy levels 
of the two branches of the one-electron band structure 
becomes very close to each other at the jFermi level. It 
has been revealed for two-Q and three-leg] Hubbard lad- 
ders that such an alignment is desirable to mimic infinite 
systems, and is indeed necessary in order to obtain corre- 
lation functions consistent with the weak-coupling theory 
at large distances. This applies to the present model as 
well. As suggested in rcf. 0, this may be an indication 
that the offset between the discrete energy levels has to 
be smaller than the spin gap in order to obtain results 
consistent with the weak-coupling theory in finite size 
systems. 

For this parameter set the Fermi wave numbers satisfy 
a relation kp — kp — 7r/2, so that the relevant pairing 
correlation becomes NNN pairs in real space according to 
eqn.y, since | exp{ikpAx) — ex'p{ikpAx)\ takes its maxi- 
mum at Ax — 2. Thus we study the pairing correlation, 



P{r) ^ (Ot^,,0,) 

Oi = (ci-fCi+2i. - Ci|Ci+2t)/\/2. 



(3) 



Figure shows the pairing correlation function for a pe- 
riodic 60-site system with U — t. The pairing correlation 
is clearly enhanced over the non-interacting case, and a 
slowly decaying component close to ~ 1/7.0.5 ^^^-^ ^^g gggj^ 
at large distances. mThe behavior is similar to that for 
two-Q and three-legQ Hubbard ladders. 

Figure is a similar plot for a larger U = 2t, about the 
largest U for which the negative-sign problem in QMC is 



overcome. In this case, the enhancement is smaller than 
that for U = t, which suggests that the exponent K is in- 
deed a decreasing function of U. However, the existence 
of the enhancement itself implies that the charge gap is 
still closed there, because otherwise the pairing correla- 
tion would decay exponentially. The exponent at large 
distances is seen to be close to unity, which indicates 
a proximity to the superconductor-insulator transition, 
so that Uc should lie somewhere above, but not too far 
from, U ^2t. 
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FIG. 1. QMC result for the pairing correlation function, 
P{r), plotted against the real space distance r in a 60-site 
half-filled Hubbard ladder with U = t a.nd t' = ~0.78t{a). 
The dashed line is the noninteracting result for the same sys- 
tem size, while the straight dashed line represents ~ l/r"'^. 
The pairing correlation function takes negative values at dis- 
tances r — S, 16, 24, which are omitted in the plot. 
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FIG. 2. A plot similar to Fig.l for U ^ 2t and t' = -0.78i. 
The straight dashed line represents ~ 1/r. 



The regime with larger U ji& exactly where we can in- 
troduce the DMRG methodO The calculation has been 
done for system sizes up to L = 28 sites with t' = — 0.8t 
here with an open boundary condition. We have kept 
up to 120 states per block with truncation errors smaller 
than 10^"* and mostly around 10^^ ^ 10^^. The charge 
gap Ac{N) and the spin gap As{N) for a half-filled A'^- 
site system are calculated by 

Ac{N) = [E{N/2 + 1, N/2 + 1) + E{N/2 - 1, N/2 - 1) 
-2£;(iV/2, iV/2)] /2, (4) 

As{N) = E{N/2 + 1, N/2 - 1) - E{N/2, N/2), (5) 

respectively, where E{N^,Ni) is the ground state energy 
for iV| up-spin and Ni down-spin electrons. We have 
checked that as far as these quantities and the parame- 
ter values adopted here are concerned, ^ 100 states per 
block suffice for the convergence. EJ 

Figure displays the charge gap for various values of 
[/ as a function of 1/L (the inverse system size). The 
results for 8- and 10-site systems are also obtained from 
the exact diagonalization in order to check the DMRG 
result. The system-size dependence is least-squares fit to 
second-order polynomials in 1/L. The extrapolation to 
L ^ 00 shows that the charge gap closes at Uc ^ 3i. 
Although the estimated Uq niay contain some errors due 
to finite-size effects, this result is consistent with the 
above QMC result that the pairing correlation exists for 
U = 2i, which serves as a lower boundary for the ex- 
istence of a spin gap. Thus the complementary QMC 
and DMRG results can indeed be combined to indicate 
a superconductor-insulator transition somewhere around 
2t<U< 3i. 
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FIG. 3. DMRG evaluation of the charge gap Ac at 
t' — —0.8t plotted against 1/L for various values of U. Dashed 
curve is a least-squares fit with second order polynomials in 
1/L. 
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FIG. 4. DMRG evaluation of the spin gap for U = St and 
t' = —0.8t plotted against 1/L. 



The existence of a metallic phase at half-filling having 
a spin gap and a dominating pairing correlation is rather 
surprising. One might be tempted to explain the ori- 
gin of the spin gap and superconductivity by regarding 
the model as essentially the same as the 't-J-J'^ model,EI 
where electrons hop with excluded double occupancies 
and interact with NN and NNN antiferromagnetic ex- 
changes. This view, however, cannot be valid, because 
we are sitting on the half-filled point, where the t-J-J' 
model is trivially insulating. In this sense the spin gap 
in the metallic phase of the present model may not be 
regarded as simply arising from a frustration in antifer- 
romagnetic exchange interactions (J, J^, which is con- 
sistent with Fabrizio's remark in ref. |lO| on the inverse 
doping effects on the spin gap between the present model 
and the t-J-J' model. 

Finally we evaluate the spin gap using DMRG. Al- 
though the spin gap is a key concept in our motivation 
mainly in the metallic phase, we present here the result 
for the insulating phase, which is of greater interest in 
the context of strongly correlated systems such as ID 
cuprates. The study of the spin gap in the metallic phase 
will be presented elsewhere. We have evaluated the mag- 
nitude of the spin gap at U = 8t with the same value 
of t' as above. In fig., we plot Ag as a function of 1/L. 
The extrapolation to the thermodynamic limit suggests 
an existence of a spin gap as large as ^ 0.16< ^ 0.17i.E3 

To summarize, a transition from a metallic state with a 
dominant pairing correlation into an insulating state has 
been shown to occur at an intermediate strength of the 
interaction of ^ 0{t). Further studies for other values of 
t' and U , and for the doped phase is under way. 
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